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ABSTRACT 

Context. The distances of pulsating stars, in particular Cepheids, are commonly measured using the parallax of pulsation technique. 
The different versions of this technique combine measurements of the linear diameter variation (from spectroscopy) and the angular 
diameter variation (from photometry or interferometry) amplitudes, to retrieve the distance in a quasi-geometrical way. However, 
the linear diameter amplitude is directly proportional to the projection factor (hereafter /7-factor), which is used to convert spectro¬ 
scopic radial velocities (i.e., disk integrated) into pulsating (i.e., photospheric) velocities. The value of the /7-factor and its possible 
dependence on the pulsation period are still widely debated. 

Aims. Our goal is to measure an observational value of the /7-factor of the type-II Cepheid k Pavonis. 

Methods. The parallax of the type-II Cepheid k Pav was measured with an accuracy of 5% using HST/FGS. We used this parallax as 
a starting point to derive the /7-factor of k Pav, using the SPIPS technique (Spectro-Photo-Interferometry of Pulsating Stars), which 
is a robust version of the parallax-of-pulsation method that employs radial velocity, interferometric and photometric data. We applied 
this technique to a combination of new VLTI/PIONIER optical interferometric angular diameters, new CORALIE and HARPS radial 
velocities, as well as multi-colour photometry and radial velocities from the literature. 

Results. We obtain a value of p = 1.26 + 0.07 for the /7-factor of k Pav. This result agrees with several of the recently derived 
Period-/7-factor relationships from the literature, as well as previous observational determinations for Cepheids. 

Conclusions. Individual estimates of the /7-factor are fundamental to calibrating the parallax of pulsation distances of Cepheids. 
Together with previous observational estimates, the projection factor we obtain points to a weak dependence of the /7-factor on period. 

Key words. Stars: individual: k Pav, Techniques: interferometric, Methods: observational, Stars: distances, Stars: variables: Cepheids 


1. Introduction 

Cepheids have been used for more than a century as standard 
candles for estimating extragalactic distances, owing to the lin¬ 
ear relationship between the logarithm of their pulsation period 
and their intrinsic magnitude: M = a (lo gP -l) + b, where P 
is the period of pulsation and M the mean absolute magnitude 
(hereafter the P-L relation). This remarkable relation was dis- 
covered empirically in 190 8 by Henrietta Leavitt ( |Leavitt|1 908; 
[Leavitt & Pickering] 19 12| and is now named after her : Leav¬ 
itt’s Law. While the slope a of the P-L relation can be estimated 
by observing a number of Cepheids located, for example, in the 
Magellanic Clouds, the zero point h is more difficult to estimate. 


* Based on observations realized with ESO facilities at Paranal Ob¬ 
servatory under program IDs 091.D-0020 and 093.D-0316. 

** Based on observations collected at ESO La Silla Observatory using 
the Coralie spectrograph mounted to the Swiss 1.2 m Euler telescope, 
under program CNTAC2014A-5. 


It was first settled by [Hertzsprung| ( ~913| ), when estimating the 
LMC distance. Thanks to the small dispersion of Leavitt’s Law, 
Cepheids have, for more than a century, been considered as one 
of the most accurate standard candles for estimating extragalac¬ 
tic distances. As a recent example, Cepheids were involved in the 
3 % m easurement of the Hubble constant Ho presented by |Riess| 
|et al.| ( [2011[ ), as calibrators of the magnitude-redshift relation of 
type la supernovae. It should be noted that half of the 3% error 
budget is attributed to the uncertainty on the Cepheid distance 
scale. The uncertain calibration of the P-L relation, as well as 
photometric biases (e.g., from reddening) can introduce impor¬ 
tant systematic uncertainties. Leavitt’s Law is also affected by 
physical effects like metallicity. A better measurement of H 0 is 
fundamental to constraining the A CDM model parameters, such 
as the number of neutrino species and the equation of state of 
the dark energy oj c . A detailed review of the methods used to de¬ 
termine the history of the expansion and their current limitations 
is given in Weinberg et al.| ( |2013| ). Better precision would also 
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help for addressing the 2<x tension between the value of |Riess| 
|et al.| |20TT ) and the value derived from Planck's CMB mod¬ 
elling ( [Planck Collaboration et al.|2014| ), and revealing potential 
biases in one of these methods. 

A better estimate of Hq requires a more accurate calibra¬ 
tion of the P-L relationship ( |Suyu et aL||2012| ). Parallax mea¬ 
surements can be used to determine the zero point of the P- 
L relationship, but they are reasonably accurate only for a few 
nearby Cepheids. Other distance estimates, such as the Baade- 
Wesselink (hereafter BW) technique (Gieren et al. [2013| ) and its 
infrared surface-bri gthness version ([Storm et al.|2011| ), the light 
echoes of RS Pup (Kervella et al. 2014), or binary Cepheid or¬ 
bital parallax ( [Pilecki et al.||2013[ |Gallenne et al.||2013| |2014| ), 
should be obtained as independent cross-checks to ensure the 
accuracy of the calibration. We are currently carrying out a long¬ 
term programme of interferometric observations of Cepheids 
in both hemispheres, using the CHARA array installed at the 
Mount Wilson Observatory in California ( [ten Brummelaar et al.| 
2005) and the Very Large Telescope Interferometer (hereafter 


VLTI) installed at the Cerro Paranal in Chile ( |Merand et al. 
|2014| ). This programme has led to the discovery of circumstel- 
lar envelopes around several Cepheids (e.g. Gallenn e~et al.|2012 


|Kervella et al.|2009~] ) and companions (e.g. |Gallenne et al. 201 3) 
The BW technique ([Baade 1926 ) combines measurements of 
the linear diameter variation (from spectroscopy) and the angular 
diameter variation (from either photometry or interferometry) to 
retrieve the distance in a quasi-geometrical way. The accuracy of 
this elegant method is although limited by the projection factor 
(hereafter /7-factor) used to convert spectroscopic radial veloc¬ 
ities into pulsating velocities representing the actual displace¬ 
ment of the photosphere. The value of the /7-factor, and its pos¬ 
sible dependence on the pulsation period, are still debated. We 
developed a dedicated software tool to estimate the distance of 
pulsating Cepheids, the SPI PS c ode ( [Merand et al.||2013| that 
we briefly present in Sect. |3.2| Taking advantage of the dis¬ 
tance to k Pavonis obtained by Benedic t et al.| ( [201 1| ) using an 
independent method, we track down in the present work the p- 
factor of this star through an inverse application of this algorithm 
(Sect. [33]). 


2. Observations and data processing 

2.1. The peculiar Cepheid k Pavonis 
As one of the closest type-II Cepheids ( [Wallerstein|2002| ), k Pav 


is classified as a member of the W Vir class ( [Rodgers & Bell 


|1963| ), which groups pulsators with periods between 10 and 20 
days. Since it is slightly brighter and bluer than the normal pop¬ 
ulation of this class, Matsunaga et al. ( 2009) propose to classify 
it as a pW st ar, which is a new category introduced by[Soszyhski 
[et al.[f2008j ) to describe those peculiar W Vir Cepheids. Type-II 
Cepheids are luminous stars, with similar behaviour to classi¬ 
cal Cepheids and variation periods in the same range. But their 
lower mass makes them representatives of the population II stars 
present in the galactic halo and the old galactic disk, k Pav also 
shows similarities with RR Lyrae, placing it at an interesting 
interface between classical Cepheids, type-II Cepheids and RR 
Lyrae pulsators. 


2.2. VLTI/PIONIER long-baseline interferometry 

Interferometric angular diameter measurements of k Pavonis 
were carried out in August and September 2013 and from June to 
August 2014 using the four-telescope beam combiner PIONIER, 
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Fig. 1 . (u,v) coverage of our PIONIER observations (each colour corre¬ 
sponds to one night of observations). The black dashed-line circle traces 
the spatial frequency where the squared visibility is equal to 0.5, assum¬ 
ing for k Pav an average diameter of 0 Avg _ - 1.182 mas (resulting from 
the present study). 

Table 1 . Main caracteristics of the calibrators used during our PIONIER 
observations. 


Star 

m v 

m H 

Oud ± cr (mas) 

HD 181019 

6.35 

3.31 

1.030 + 0.014 

HD 171042 

7.51 

3.63 

1.030 + 0.014 

HD 175782 

6.80 

3.05 

1.337 + 0.019 

HD 186530 

7.41 

4.01 

0.896 + 0.012 


Notes. Calibrators were selected from the catalogue of [Merand et al.| 
( [2005a[ ). We indicate the uniform disk diameter in H band and the cor¬ 
responding uncertainty. 


installed at the VLT Interferometer in northern Chile (Berger 


et al.|2010t|Le Bouquin et al.|2011| ). The observations were un¬ 


dertaken in three spectral channels of the H band (1.59 pm, 1.67 
pm, and 1.76 pm), corresponding to a low spectral resolution 
of R = 40. We used the four relocatable 1.8-metre Auxiliary 
Telescopes (ATs), installed at stations A1-G1-J3-K(Q This con¬ 
figuration offers the longest available baselines (from 57 meters 
between the stations K0 and J3, up to 140 metres between Al 
and J3) that are required to get the necessary angular resolution 
to resolve the apparent disk of k Pav (6 ~ 1.2 mas). 

Figure [I] shows a map of the ( u, v) plane coverage of our ob¬ 
servations. We obtained a very good sampling of the pulsation 
curve, including data points close to the maximum and the min¬ 
imum diameters. That is particularly interesting for the present 
work, because the most precise diameter variation is needed to 
derive an accurate value of the /7-factor. The most suitable cal¬ 


ibrators have been chosen from the catalogue of [Mer and et al. 
( 2005a| ). Their main characteristics are given in Table [I] 


1 https ://www.eso.org/paranal/telescopes/vlti/configuration/ 

















































































J. Breitfelder et al.: Observational calibration of the projection factor of Cepheids 



Fig. 2. Squared visibilities measured with PIONIER at the minimum 
and the maximum diameters. The data are fitted with a uniform disk 
model leading to the diameters of 0 min = 0.976 mas (MJD=56 539, in 
yellow) and 6 max = 1.273 mas (MJD=56 871, in green). Errorbars have 
been removed for clarity. 


Table 2. PIONIER observations of k Pav. We give here the mean MJD 
(defined by JD-2 400 000.5) of each observing night, the correspond¬ 
ing phase (taking 0 O at the maximum of luminosity), the best uniform 
disk diameter adjusted on the squared visibility measurements and its 
uncertainty. We also give the x 1 resulting from the fit of the squared 
visibilities with a uniform disk model. 


MJD 

Phase 

@UD — t^stat. — ^syst. (mas) 

* 2 

56508.024 

0.375 

1.256 ±0.004 ±0.014 

0.96 

56509.009 

0.483 

1.243 ±0.002 ±0.015 

3.30 

56538.065 

0.682 

1.101 ±0.012 + 0.016 

1.30 

56539.046 

0.790 

0.976 ± 0.016 ± 0.016 

1.35 

56540.020 

0.898 

1.159 ±0.012 ±0.015 

0.59 

56831.126 

0.948 

1.101 ±0.004 + 0.014 

2.21 

56833.166 

0.173 

1.262 ±0.002 ±0.014 

1.73 

56866.992 

0.897 

1.057 ±0.002 ±0.014 

1.01 

56869.122 

0.132 

1.222 ±0.007 ±0.014 

1.48 

56871.055 

0.344 

1.273 ±0.004 ±0.014 

0.94 

56874.090 

0.679 

1.111 ±0.004 ±0.014 

1.10 

56894.017 

0.872 

1.026 ±0.001 ±0.014 

0.69 


The raw data were reduced using the pndrs data reduction 
software of PIONIER (Le Bouquin et al. 2011| , which produces 
calibrated squared visibilities and phase closures. The resulting 
01 FITS format of the data allowed us to compute a uniform disk 
(UD) value for each observing epoch wit h the model fitting soft¬ 
ware LITPrc[^ ( Tallon-Bosc et al.| |2008). The results are listed 
in Table [2] For each night of observations we give the mean Ju¬ 
lian date, the phase, the UD diameter, and its uncertainty. The 
error given for each diameter includes a statistical error given by 
the model fitting software and a systematic one (which actually 
dominates the error bugdet), because of the uncertainty on the 
calibrator diameters. We therefore define this systematic error as 
being the mean of the errors of all the calibrators. The phases 
given in Table [2] were computed by using the ephemeris derived 
from the O-C diagram shown in Fig. [3] (see details in Sect. |2.3| ). 
The variations in the diameter of the Cepheid can be appreciated 
in Fig. [2} which represents the PIONIER squared visibilities and 
the corresponding best UD model for the maximum and mini¬ 
mum diameters. 


2.3. Period changes, overall phasing, and photometry 

The pulsation period of k Pav shows large and fast variations 
that complicate the phasing of datasets from different epochs. 
By applying the statistical Eddington-Plakidis method (Edding- 
ton & Plakidi s~p929| ) on a large photometric dataset, |Berdnikov 
& Stevens| (|2009 ) showed that the period variations of k Pav 
are erratic, which has already been suggested in previous stud¬ 
ies (e.g. |Feast et aLl|2008[ |Wallerstein et al.|[l992| ), with a rel¬ 
atively high degree of randomness. To help the phasing of the 
data used in the present study, we computed ephemerides from 
the O-C diagram shown in Fig. [3] When constructing the O- 
C diagram, only photoelectric and CCD photometric data have 
been taken into account. Depending on the number of observa¬ 
tions and phase coverage of the individual datasets, a weight of 
1, 2, or 3 has been assigned to the derived moment of the nor¬ 
mal maximum brightness. In Fig. [3] the size of the filled circles 
refers to the weight assigned to the given residual (O-C differ¬ 

2 The LITpro software is available at http://www.jmmc.fr/litpro 


kappa Pavonis C = 2450373.2847 + 9.087295E 



JD 2400000+ 


Fig. 3. O-C diagram of k Pav. A weighted least squares fit to the resid¬ 
uals leads to T 0 = 2450373.2847 and P = 9.0873. Considering only the 
data with JD>2450000, we get T 0 = 2450374.2938 and P = 9.0827. 


ence). A weighted least squares fit of the data for JD>2440000 
leads to To = 2450373.2847 and P = 9.0873 days. Considering 
only the data for JD>2450000, we get To = 2450374.2938 and 
P = 9.0827 days. These values have been used to phase our most 
recent data. 


We selected the following photometric data from the liter¬ 
ature: Hipparcos and Tycho photometry from the |ESA| ([1997) 
catalogu e (see also|van Leeuwen et al.| ( |1997| )), JHK photome- 


try from ] 

Feast et al. 

(2008) (hereafter F08), VBLUW photome- 

try from 1 

^Valraven et al. (1964), and UBVR C I C photometry from 

Berdniko 

v (2008). We divided this last dataset into three groups 


covering different epochs of around 700 days and separated by 
around 350 days (Group 1: MJD from 50347 to 50917, Group 2: 
51248 to 51972 and Group 3: 52323 to 53118). We then phased 
the resulting photometric sequences separately. Considering the 
erratic changes in the period of k Pav, we computed specific 
(P, To) ephemeris for each dataset. We used the SPIPS code to 
derive a reference epoch and the corresponding pulsation period. 
In the case of the Hipparcos and Tycho data, we introduced a 
linear variation in the period, which allowed us to reach a more 
satisfying phasing. For each photometric dataset, the reference 
epoch was set at the maximum luminosity of the star, and taken 
as close as possible to the centre of the epoch range. The result- 
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Table 4. Radial velocity measurements deduced from the spectra ob¬ 
tained with CORALIE and HARPS, both installed in La Silla obser¬ 
vatory in Chile. We adopted a total error budget of 500 m/s for each 
observation. 


MJD 

Phase 

Vrad ± O'tot. (km/s) 

Instrument 

56750.416 

0.096 

28.74 ± 0.50 

CORALIE 

56751.412 

0.206 

33.32 ± 0.50 

CORALIE 

56826.428 

0.465 

44.54 ± 0.50 

CORALIE 

56827.425 

0.575 

49.31+0.50 

CORALIE 

56876.174 

0.942 

26.09 ± 0.50 

HARPS 

56876.976 

0.031 

27.55 ± 0.50 

HARPS 

56877.974 

0.140 

31.77 + 0.50 

HARPS 

56878.975 

0.251 

36.26 + 0.50 

HARPS 

56907.984 

0.444 

44.13 + 0.50 

HARPS 

56908.983 

0.554 

48.92 + 0.50 

HARPS 

56909.994 

0.666 

51.72 + 0.50 

HARPS 


ing periods and reference epochs MJD 0 are given in Tableland 
were used to compute the phases of the observations. 

We note that the Hipparcos and Tycho data seem to be rela¬ 
tively dispersed around the phase 0.6, given their low errorbars. 
This could be due to possible amplitude variations, which have 
already been observed in the photometry of type II Cepheids and 
overtone pulsators ([ Evans et a l.j| 2014| . However, this does not af¬ 
fect the results of the present study. 


Table [3] also sums up the method used to phase the other ob¬ 
servables : for the ra dial v elocities from [Waller stein etaL] ( |1992| ) 
(see details in Sect. |2.4| ), we kept the phases provided by the 
author and shifted the whole curve to make it match the photom¬ 
etry. Our CORALIE and HARPS data were phased by using the 
ephemeris derived from the O-C diagram (see Fig. [3}, and then 
shifted into phase agreement with the rest of the data. We used 
the same ephemeris for our PIONIER diameters, but we did not 
need to add any phase shift. 


2.4. Radial velocities 


We retrieved three sets of radial velocity measurements from the 
literature, resulting from o bservatio ns carried out betwee n 1904 
and 1 918 (|Jacobsen|1929|), in 1961 ( [Rodgers & Bell|1963] ) and in 
1988 ([Wallerstein et al.|1992|) (hereafter W92 ). As the data from 
| Jacobsen ( [1929 ) and Rodg ers & Bell| ([1963 ) show a relatively 
high dispersion, we only used the metallic-line radial velocities 
from W92. The convention used to derive the phases given in 
W92 is uncertain (Wallerstein, private communication) and by 
computing our own ephemerides, we did not succeed in obtain¬ 
ing a radial velocity curve that was as well-phased as the one 
given in W92. We therefore decided to keep the phases given 
in the paper. Because the zero-phase definition in W92 does not 
correspond to the maximum flux, we shifted the original phase 
values to match the zero-phase convention of the photometry. 

We also obtained contemporaneous and very precise radial 
velocities (RVs) between November 2013 and June 2014 using 
the high-resolution Echelle spectrographs HARPS mounted on 
the ESO 3.6m telescope and CORALIE mounted on the Swiss 
1.2m Euler telescope, both of which are located at ESO’s La 
Silla observatory in Chile. The RVs were derived using the cross¬ 
correlation technique (Baranne et al.|1996[|Pepe et al.|2 002) and 
a standard Gaussian fit. 
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Unfortunately, the new measurements are not sufficient to 
determine a very precise pulsation period. Therefore, we com¬ 
puted the phases of the RV measurements using the ephemeris 
used for the PIONIER data. We then shifted the whole curve by 
adding a constant in phase, until obtaining the lowest reduced 
X 1 at the end of the fitting process. The phasing process is sum¬ 
marized in Table [3] While we took great care to obtain the best 
possible phasing for the new RV data, there is a remaining un¬ 
certainty regarding the phases of the new measurements. This 
problem is exacerbated by know n random period fluctuations 
(see details given in section |2.3| ). To reduce the sensitivity of 
our method to phase errors of the new RV data (this easily incurs 
errors of several 100 m/s and thus dominates the uncertainty of 
the RV curve), we attribute a reduced weight to the new RV data 
by adopting a fixed error budget of 500 m/s in the fit. This does 
not constitute a limiting factor for our results, since the accuracy 
of our p-factor determination is limited by the parallax accuracy 
of 5% (see details given in section [33] ). 


3. The SPIPS algorithm 

3.1. The parallax of pulsation method and the p-factor 
limitation 

The fundamental equation of the parallax of pulsation method 
can be written as follows, where Q t= o is the angular diameter at 
the maximum of luminosity, d the distance and v pu i s the pulsation 
velocity of the atmosphere: 

0(t) = 0 t =o + 2 Vpuis ( t)c/t ( 1 ) 

The main limitation of the parallax of pulsation technique 
comes from the projection factor p used to convert the radial ve¬ 
locities from the Doppler shift of spectral lines (a disk-integrated 
quantity) into pulsation velocities (the displacement velocity of 
the photosphere over the pulsation cycle): 


^puls — P Uad 


( 2 ) 


We do not consider here the effects of amplitude modulations 
as recently reported for Cepheids by lAnderson^O^) . We can 
therefore rewrite the main equation of the parallax of pulsation 
as follows: 


0(t) 


d 


Ro+ f 

, Jo 


PVmdWdrj 


(3) 


Only the projected component of the velocity on the line of 
sight contributes to the measured Doppler shift of the spectral 
lines. The spherical geometry of the star results in a projection 
effect corresponding to a value of p - 1.5 for a uniform bright¬ 
ness sphere. The limb darkening and the dynamical behaviour of 
the line-forming regions in the stellar atmosphere are expected 
to reduce the value of p below this value ( [Nardetto et al.|2007[ ). 
Unfortunately, these effects can hardly be quantified separately. 

The combination of these different effects results in a rela¬ 
tively large dispersion of the p-factor estimates found in the liter¬ 
ature. Some recent studies propose either a constant value ( Feast 
|et al.|2Q08]|Groenewegen|20Q7]|Merand et al.|2005b|) or Period 


p-factor relationshi ps ([Gro enewegen 2013; Stor m et aL||2012[ 

Nge ow et al.||2012| Nardetto et al.||2009| ). As a consequence of 
this uncertainty, the truly unbiased quantity that can be derived 
using the parallax of pulsation method is the ratio d/p, where 
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Table 3. Ephemerides used to phase the interferometric, spectroscopic, and photometric data. 


Data reference 

Phase shift MJD 0 

Period (days) Period drift (s/year) 



Radial velocities 


Waller stein et al. 

<N 

OS 

OS 

) 0 + 0.171 

- 

Present work 

0 + 0.034 50373.793 

9.0827 



Interferometry 


Present work 

50373.793 

9.0827 


Photometry 


Berdnikov ( 

2008 

(i) 


- 

50646.541 

9.0897 

- 

Berdnikov ( 

2008 

(2) 


- 

51609.425 

9.0804 

- 

Berdnikov 

2008 

(3) 


- 

52353.321 

9.0669 

- 

Feast et al 

.(2008 

) 


- 

46800.743 

9.0789 

- 

Walraven et 

al. ( 

19( 

)4 

) 

- 

37913.807 

9.0833 

- 

ESAU 

997 

') 




- 

47482.679 

9.0888 

65.895 


d is the distance (Eq.[3]). Observational determinations of p are 
therefore critical for better constraining the p-factor models and 
eventually establishing the parallax of pulsation method on solid 
foundations. A review of the role and importance of the ^-factor 
can be found in [Nardetto et al.[ ( [2014| ). 

The observational calibration of the ^-factor requires in¬ 
dependent measurements of the distance d for a sample of 
Cepheids, to waive the degeneracy of the d/p ratio produced by 
the parallax of pulsation method. In this case, when the distance 
is known (e.g. from trigonometric parallax, light echoes, or bi¬ 
nary Cepheids), the inversion of the method gives access to the 
p-factor. T his has already been done on the classical Cepheid 5 
Cephei by |Merand et al.| ( [2b05b| ), yielding p = 1.27 ± 0.06. We 
apply here this method to k Pav using the enhanced SPIPS al¬ 
gorithm (Sect. [3]) and considering the accurate parallax recently 
measured by [Benedict et al.| ( [2011| ) using the HST/FGS interfer¬ 
ometer (7r = 5.57 ± 0.28 mas). 


3.2. Overview of the algorithm 


Classical BW methods are limited by various systematic er¬ 
rors (e.g. photometry biases or reddening) that affect either 
the photometry or the spectroscopic observables. To overcome 
this, we developed a dedicated tool (SPIPS for Spectro-Photo- 
Interferometry of Pulsating Stars, Mera nd et al.| ( [2013] ); |Breit- 


feld er et ah] ( |2014| ); Merand et al., in prep.) that computes a 
global fit of all the available data (i.e., radial velocities, inter¬ 
ferometric squared visibilities, spectroscopically determined ef¬ 
fective temperatures (T e ff), colours and magnitudes in various 
bands and filters). This combination of different observables al¬ 
lows us to reach a better accuracy on our measurements (2 to 5 % 
uncertainty on the distance for an individual Cepheid). The par¬ 
tial redundancy of the data (e.g., interferometry and atmospheric 
models applied to the photometry to estimate the angular diam¬ 
eter) results in much improved robustness of the fitting process. 
Besides that, the integration of physical models in our code (e.g. 
ATLAS9 atmospherical models from lKurucz| ( [l979] )) reduces the 
statistical biases (owing to the calibration of the zero point of 
the photometric filters, for instance). The larger overall amount 
of observational data also reduces the statistical errors on the re¬ 
sulting parameters (e.g., the distance, the colour excess E(B-V ), 
and the mean angular diameter and effective temperature). The 
X 1 minimization process is optimized by defining the global x 2 
as the average of the specific reduced x 1 values for each data 



pulsation phase 

Fig. 4. SPIPS code applied on the type-I I Cepheid k Pav. Above: Radial 
velocities from |Wallerstein et al.|jl992| (the typical error bar is shown 
in the bottom left corner) and new measurement from CORALIE and 
HARPS (the size of the points corresponds to the errorbar), fitted using 
spline functions. Below: Uniform disk angular diameters deduced from 
our PIONIER interferometric observations. For each panel we indicate 
the reference and the individual reduced^ 2 . 


set. This allows us to adjust the relative weight of each type of 
data in the global fit. Otherwise, in the present case, the diame¬ 
ter and the radial velocity would not contribute as much as the 
photometry in the minimization process. 
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Fig. 6. SPIPS code applied on the type-II Cepheid k Pav : magnitudes and colours fitted with Fourier series. For each dataset we precise the filter, 
the reference, and the individual reduced^ 2 . 



Fig. 5. SPIPS code applied on the type-II Cepheid k Pav. Effective tem¬ 
perature curve deduced from ATLAS9 atmospheric model grids (|Castelli| 
|& Kurucz|2004l |Kurucz|2005] > 


3.3. Projection factor 

We fitted the radial velocity with spline functions defined by four 
adjustable floating nodes. Although it is numerically less stable 
than Fourier series, it is necessary to reproduce the strong slope 
of the radial velocity curve that occurs between the phases 0.8 
and 1.0 without introducing high frequency variations in the rest 
of the curve. The photometry curve is adjusted to the data using 
a second-order Fourier series. We shifted the (B - R ) and (B -1) 
colours from [Berdnikov] p008| ) vertically, by subtracting 0.008 
magnitudes to (B - I) and 0.038 magnitudes to (B - R). Other¬ 
wise, an offset always remained between the data and the model, 
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probably because of a bias introduced by the calibration of the 
zero point of the filter used for the observations in I and R. It is 
then important to note that only the shape of these two particular 
curves continues to be constraining. To give an equivalent weight 
to the different observables in the fitting process, we multiplied 
the errors by the normalization factors (NF) given in Table [6] 
These coefficients depend on the number of points (N po i nts ) in 
each dataset. Table [b] also gives the individual reduced^ 2 of the 
different datasets. 

We fixed the metallicity at [Fe/ H] = 0.0 ([Luck & Bond 


1989| and the distance at d = 179 pc (Benedict et al .|2011| ). The 


(Be 

best-fit parameters are given in Table [5] The values and statisti¬ 
cal uncertainties have been determined through a Monte-Carlo 
procedure. We obtain for k Pav a /7-factor p - 1.26 ± 0.04 stat ± 
0.06 syst , with a systematic error given by the limited 5% preci¬ 
sion on the distance. The final value is therefore p - 1.26 ± 0.07. 
The systematic error on the average T q q has been deduced from 
the Stefan-Boltzmann law, after considering a systematic error 
of 5% on the photometry and of 1.4% on the angular diameter 
(given by the limited precision of the calibrators diameters). The 
systematic uncertainty on the reddening has been estimated by 
computing the maximum and minimum values of the colour in¬ 
dex (B - V)o obtained in the uncertainty range of the effective 
temperature ( Flower|1996[ ). 


4. Discussion 

The ^-factor value of k Pav that we obtain ( p = 1.26 ± 0.07) 
is significantly higher than the value of 0.93 ±0.11 proposed 
by F08. A /7-factor smaller than unity would imply that the 
limb darkening of the star is extremely high, and would gen¬ 
erally not have a very clear physical explanation. The present 
value, however, agrees well with recently published Period-/7- 
factor relationships, which give for k Pav (P = 9.09 d) /7-factors 
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Table 5. Best-fit output parameters given by the SPIPS code. 


Parameter Value <x stat . <x sys . 


^-factor 

0 (UD) at 0 = 0 (mas) 
Avg. 6 (UD) (mas) 
v r (km/s) 

E(B - V) (mag) 

Avg. r eff (K) 

Avg. radius ( R e ) 

Final reduced^ 2 


1.26 

0.04 

0.06 

1.1654 

0.0025 

0.014 

1.1823 

0.0021 

0.014 

37.87 

0.18 

0.50 

0.02 

0.01 

0.04 

5739 

9 

107 

22.83 

2.62 

0.04 

1.14 


Table 6 . Weighting of the different datasets and individual reduced^ 2 . 


Observable 

Ref. 

Npoints 

NF 

+ 

V Johnson 

(a) 

162 

1.546 

1.42 

U Walraven 

(b) 

16 

0.486 

1.43 

V Tycho 

(c, d) 

163 

1.550 

4.37 

B Tycho 

(c, d) 

163 

1.550 

6.45 

HP £2000 

(c, d) 

182 

4.914* 

6.97 

JCTIO 

(e) 

25 

0.607 

1.46 

Hctio 

(e) 

25 

0.607 

0.51 

Kctio 

(e) 

25 

0.607 

1.94 

^Walraven ~ Bw air avert 

(b) 

16 

0.486 

8.08 

B Johnson ~ V Johnson 

(a) 

56 

0.909 

1.51 

V Johnson ~ IBessel 

(a) 

164 

1.555 

0.57 

V Johnson ~ RBessel 

(a) 

27 

0.631 

0.95 


CD 1 

4 

0.243 


Vrad 

CO 2 

7 

0.321 

1.43 


(g) 

38 

0.749 


6 (UD) 

(f) 

12 

0.421 

2.30 


Notes. References : (a) Berdnikov (2008); (b)|Walraven et al. 

(1964); 

(c) ESA (1997); (d) van Leeuwen et al. (1997); (e)|Feast et al. 

(2008); 

(f) Present work ^CORALIE and Z HARPS); (g) Walraven et al. 

(1964). 

*The normalisation factor of the Hipparcos data has been multiplied by 


three to balance the very low uncertainties of the measurements, which 
would otherwise give to this dataset too much weight in the fitting pro¬ 
cess. 


of 1.27 ( |Groenewegen||2013|), 1.29 + 0.06 ( [Ngeow et al.|2012| ), 
and 1.23 + 0.10 ( [Nardetto et al.|20Q9| ). It i s also consistent with 
the empirical measurements obtained by |Pilecki et al.| (|2013|) 
on the LMC Cepheid OGLE-LMC-CEP-0227 (P = 3.90d, 
p = 1.21 ± 0.05), and by [Merand et al.|| ( |2005b| ) on £Cep 
(P = 5.37 d, p = 1.27+0.06). Our /7-factor, however, diffe rs from 
the p = 1.37 an d p = 1.359 + 0.003 values predicted by [Storm] 


let al. 


( |2012| ) and Nei lson et al.| ( |2012| ) respectively. The average 


angular diameter is in good agreement with the value derived 
by Gallenne et al.| ( |20T2] ) (#ud = 1.04 + 0.04 mas at 0 = 0.9). 
Converted into linear radius for a distance of d = 179 pc, we 
obtain an average photospheric radius of 22.8 R 0 . |Balog et ak| 
( 1997| suggest a comparable value of 19 + 5R 0 , also derived 
using the parallax of pulsation method. The average SPIPS ef¬ 
fective temperature of T q q = 5739 K is slightly higher than the 
typical values found in the literature. In particular, [Luck & Bond] 
(1989) find5500K, and |Gallenne etaL] ( [2QT2l ) find T eff = 5750K 


at 0 = 0.94, which corresponds to ~ 6336 K at the same phase 
in the present study. We find an extinction comparable to the 


value suggested by F08 ( E(B - V) = 0.017). It is important to 
stress that the relevance of this parameter relies on the choice 
of a reddening law, which is in the present case the reddening 
law from |Fitzpatrick[ ( j!999j ), taken forR v = 3.1 and differs from 
the methods used in F08. However, the systematic errors of both 
measurements dominate these low extinction values. 

It was suggested that k Pav belongs to a binary system ( [Feast 
et al.|2008] ). The contribution of a stellar companion could have 
a non-negligible influence on the photometry and radial velocity 
of the star. We checked our interferometric dataset for the possi¬ 
ble presence of a secondary component, by fitting a binary star 
model that takes the PIONIER closure phases into account. This 
fitting technique is very sensitive to the presence of companions, 
with a contrast exceeding 100:1 or more (Ab sil et al.|2011||Gal-| 
lenn e et al.|2013||2014| ), but we did not identify any companion 
of k Pav. Our data allowed us to obtain an upper limit of 1% at 
5<x on the flux ratio between the two components. This result has 
been derived from the longest observing sequence (about 3 hours 
of science and calibrators alternations), which allows reaching a 
good sensitivity and better uv coverage. Considering the low flux 
ratio limit, it is unlikely that the data used in the present study 
could have been biased by the presence of a companion. 

The SPIPS code also allows us to consider a possible in¬ 
frared excess in the fitting process, to track the possible pres¬ 
ence of a circumstellar envelope. In the case of k Pav we find an 
excess of 4.5 + 0.5 % in the K band. However, it does not im¬ 
prove the quality of the fitting process significantly, so we prefer 
not to conclude anything about the presence of an actual excess. 
We nevertheless underline that a circumstellar envelope has been 
found by Gallenn e~et al.| ( |20T2] ), who identified an infrared ex¬ 
cess of about 20% between lO/xm and 20 pm. 


5. Conclusion 


We report the first observational measurement of the projection 
factor of the type-II Cepheid k Pav. We combined the HST/FGS 
parallax from Benedi ct et al.| ( |2011| ) with new interferometric 
observations from the VLTI/PIONIER instrument, and an ex¬ 
tensive set of radial velocities and photometry. Because the pe¬ 
riod of the star shows unpredictable variations on relatively short 
timescales, a careful phase adjustment was required to phase 
the different observing epochs properly. We obtained a value of 
p = 1.26 + 0.07, which agrees with the Period-/7-factor rela¬ 
tion proposed by |Nardetto et al.|(j2009| and with the empirical 
measurement obtained by Pilecki et al. (2013} in the LMC. It 
is also consistent with the /7-factor measured by |Merand et al. ' 


( j2005b[ ) for S Cep. Although the range of periods presently cov¬ 
ered by observational /7-factor measurements is still limited, our 
result points at a relatively weak dependence of the /7-factor on 
the period, because short- and intermediate-period type-II and 
classical Cepheids likely share the same /7-factor within « 5%. 
Observational measurements of the /7-factor are difficult, but also 
essential for the calibration of the distance scale. This factor is 
presently one of the most important fundamental limitations on 
the accuracy of the parallax of pulsation distances used to cal¬ 
ibrate period-luminosity relationships of Cepheids. Radial ve¬ 
locity modulations (though not evident in the present case) can 
also lead to systematic errors, as discovered recently by [Ander¬ 


son] ( |2014| ). Observational estimates of this parameter are also 
essential for constraining the /7-factor models. The Gaia satellite 
is currently measuring accurate parallaxes for a large number 
of Galactic Cepheids. This will enable a thorough study of the 
dependence of the /7-factor with period and other stellar param- 
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eters, and provide us with a solid, unbiased calibration of the 
parallax of pulsation technique. 
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